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ABSTRACT 

Eight years of optical photometry from OGLE-III are presented of the optical 
counterpart to the High Mass X-ray Binary system, SMC X-l. The optical data pro- 
vide the best view to date of the modulation of the light from the system at the binary 
period of 3.9d. In addition, it is shown for the first time that the light is also modu- 
lated at the superorbital period - a period associated with the precessing of a warped 
accretion disk around the neutron star partner. The implications for the sources of 
the various components of the optical light in this system are discussed. 
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1 INTRODUCTION AND BACKGROUND 

A survey of the literature reveals ~200 High Mass X-ray 
Binaries (HMXBs) known in our Galaxy and the Magel- 
lanic Clouds (Liu et al., 2005, 2006). Of these, there is a 
substantial population of HMXBs in the Small Magellanic 
Cloud (SMC) comparable in number to the Galactic popula- 
tion (Coe et al., 2009). Of the spectrally classified HMXBs 
in the Galaxy there are 21 supergiant systems and 28 Be 
star systems in the current literature. However, unlike that 
Galactic population, all except one of the SMC HMXBs are 
Be star systems. The one exception is the supergiant HMXB 
system SMC X-l which is known to be unusual in having a 
rather low mass neutron star in the binary partnership (Val 
Baker, Norton, & Quaintrell 2005). 

SMC X-l is a binary system consisting of a neutron star, 
with pulse period 0.71s (Leong et al. 1971), and the BOI 
super-giant companion Skl60 (Webster et al. 1972; Liller 
1973). The system has an orbital period of 3.89d and a vari- 
able superorbital period, (ranging from 40 to 70 days) with 
an average value of 55. 8d (Trowbridge et al. 2007). A warped 
accretion disk which partly obscures the neutron star when 
periodically passing in front of it is thought to be responsi- 
ble for the superorbital modulations (Katz, 1973), and the 
variation in the superorbital period is explained by Clarkson 
et al. (2003) as due to multiple warp modes in the disk. 

In this paper the relationship between the X-ray mod- 
ulations and the optical light is explored by making use of 
high-quality OGLE-III photometric monitoring over eight 
years. Detailed modelling of the binary modulation permits 
explicit confirmation of the unusually low mass of the neu- 



tron star. In addition, for the first time, evidence is presented 
for the presence of an optical modulation at the super-orbital 
period of the system. 



2 OPTICAL DATA 

Optical photometric data of the SMC X-l system were ob- 
tained from the OGLE-III project (Udalski, 2003). The I- 
band data were collected approximately every night during 
the SMC observing season for eight years (April 2001 to 
June 2009). These data are illustrated in Figure [1] It is im- 
mediately apparent from inspecting this figure that there is 
an underlying variability of the data with a range of ±0.1 
magnitudes. 



3 X-RAY DATA 

The ASM website presents the X-ray data in two different 
formats: daily averages - reliable data with small errors; and 
'dwell by dwell' - a 90s exposure data point for every ninety 
minute long orbit which, individually have larger errors. Af- 
ter experimenting with both formats the 'dwell by dwell' 
data proved more appropriate for this work due to the in- 
creased sample size. These data are shown in Figure [2] 
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Figure X. OGLE-III I- band photometry of the optical counter- 
part within the SMC X-l system. 
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Figure 2. RXTE/ASM dwell-by-dwell X-ray data for SMC X-l. 
The source is clearly detected for over 13 years. 

4 DATA ANALYSIS 
4.1 The orbital period 

Lomb-Scargle periodograms (Lomb 1976; Scargle 1982) were 
produced from the ASM and OGLE data. The orbital pe- 
riod is a prominent peak in both the X-ray and optical power 
spectra, see Figure [3] • The second highest peak in the op- 
tical data is the product of the window function - a pe- 
riodicity introduced by the daily sampling method of the 
OGLE data. The binary period found in the ASM data was 
3.89229090±0.0000003 days - in good agreement with the 
value presented in Wojdowski et al. (1998) - see below. 
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Figure 3. Lomb-Scargle power spectra for SMC X-l over the 
period range 1.25 — lOd. The power spectrum from the X-ray 
data is shown in the upper panel and that from the optical data 
in lower panel. The orbital period and its harmonics are indicated 
in the top panel. 



The next step was to fold both optical and X-ray data 
sets on the orbital period. However in the case of SMC X- 
1 the orbital period is not constant: it is decaying (Woj- 
dowski et al. 1998) and the time of central eclipse drifts over 
the years. The centre time of the Nth eclipse as defined by 
Wodjowski et al. is given by tN — a + a.iN + C12N 2 , where 
a defines the initial zero-phase point: 2442836.18278(20) 
(days), ai is the orbital period: 3.89229090(43) (days), a 2 = 
—6.953(28) x 10 -8 and N is the number of orbits since ao- 
The additional 'squared' term compensates for the orbital 
decay. A barycenter correction was also applied to the data 
and the resulting folded X-ray and optical (I band) light 
curves are shown in Figure [4] 

The neutron star eclipse is the dominating feature in 
the upper panel of Figure H the RXTE/ASM folded light 
curve, manifesting between phases 0.85 and 1.15. 

Some considerable effort was made to reproduce the de- 
tails shown in the equivalent figure (Fig. 1) in Trowbridge 
et al (2007). Data from the same time period were used and 
data points were removed for which the RXTE/ASM solu- 
tion gave x 2 > 1-5 (as was done by Trowbridge et al.), and 
then values > 1.3. Various combinations of the RXTE/ASM 
'A', 'B' and 'C band data (C is the hardest band, A the soft- 
est) were tried, but without success. Unfortunately, we failed 
to replicate their light curve and, in particular, were unable 
to see the mid-eclipse X-ray 'bounce-back' feature that they 
report in their work. 

It was noted that if we included all RXTE/ASM data 
in our light curve then the scatter on each point was greater 
than if only data from the time period reported in Trow- 
bridge et al (2007) were used. This is probably due to the 
gradual deterioration of the RXTE/ASM over the years. 
Therefore only X-ray data from this time period were used 
in this work. 
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Figure 4. Data folded at the binary period - see text for details of 
the binary period. X-ray data (upper panel) from RXTE/ASM in 
units of ASM counts/s. Optical data (lower panel) from OGLE- 
III is in flux units of 10 — 15 erg cm~ 2 s — 1 A" 1 and folded with the 
same orbital parameters as the X-ray data. 



Since the X-ray flux of SMC X-I varies by almost 100% 
over the course of each superorbital cycle, in an attempt to 
improve the quality of our folded light curve only data near 
the times of maximum flux were used i.e. when the superor- 
bital modulations were at the peak of their cycles. Upon pro- 
ducing this new folded light curve, however, the base level 
sections of the light curve were improved but the eclipse 
suffered from the reduced sample size (which was now three 
times smaller than it was originally). Other combinations of 
superorbitally selected data were selected in an attempt to 
further refine the folded light curve, but ultimately decided 
against - not wanting to over-manipulate the data. 

The optical light curve shown in the lower panel of Fig- 
ure |4] exhibits the classical ellipsoidal variations of a binary 
system. What is not clear, however, is the asymmetrical na- 
ture of the light curve: one would expect modulations to be 
equal, peaks and dips to reach the same maximum and min- 
imum flux value. This is has been interpreted successfully 
by Rawls et al (2011) in terms of X-ray heating, ellipsoidal 
variations and a possible accretion disk. This fit is revisited 
and addressed further in the Discussion section of the paper. 

Following on from the shape of the folded I band light 
curve, it was decided to include the less-sampled V band 
data from OGLE-III, with the intention of creating a folded 
colour (V-I) phase diagram. In an attempt to maximise the 
signal-to-noise ratio the possibility was explored of includ- 
ing more V band data from Penfold and Warren (1975) and 
Paradijs (1977). However, these extra data were eventually 
discarded as they only degraded the much superior quality 
of the OGLE-III data. This meant sacrificing sample size, 
hence a lower phase resolution colour light curve of only fif- 
teen phase bins. This was produced and this is shown in 
Figure [5] This colour graph reveals an increase in temper- 
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Figure 5. Ratio of the flux in the OGLE-III V-band to that in 
the I-band as a function of binary phase. 

ature at the time of eclipse (phase 0.0) and a decrease in 
temperature as the X-ray source passes in front of the star. 



4.2 The superorbital period 

Once again Lomb-Scargle periodograms were produced from 
the RXTE/ASM and OGLE-III I-band data, this time cov- 
ering the period range reported for the superorbital period 
- see Figure [6] In the X-ray data a strong but complicated 
peak is revealed at 55.8 days (the average superorbital pe- 
riod) with the power of ~200. This same period is also 
present in the OGLE data, albeit with something of a lower 
Lomb-Scargle power - a modest ~10. 

A larger peak at 50. 7d can be seen to the right of the 
optical superorbital spike in the lower panel of Figure]!)] The 
origin of this feature is unclear. A randomising test was done 
on the optical data which returned L-S power values of 12.4, 
15.0 and 17.4 at 99%, 99.9% and 99.99% levels of confidence 
respectively. So this peak is significant at the ~ 99.9% level. 
One possibility is stellar non-radial pulsation with a period 
of 0.978 days beating with the one day sampling to pro- 
duce this feature. There are several examples of this effect 
reported in Bird et al.(2012) in their study of the OGLE 
data of another similar class of HMXBs - Be/X-ray binary 
systems. 

Nonetheless, the possibility that this period could be re- 
lated to the superorbital period was explored by splitting the 
OGLE data into three equal length sections and producing 
periodograms of each. If this period wandered around fol- 
lowing the original superorbital period, then Lomb-Scargle 
of each third should look different. While the true superor- 
bital period does indeed vary, this new peak remains at the 
same approximate period throughout - offering support to 
the non-radial pulsation model. 

Despite the low power of the secondary peak at 55. 8d, 
however, we can be confident that some of the optical light 
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Figure 6. Lomb-Scargle periodograms of the X-ray (upper panel) 
and optical (lower panel) data sets showing the region around the 
super-orbital period. The vertical dotted line in the lower panel 
indicates the peak position from the upper panel and represents 
a period of 55. 8d. 



is being modulated on the same superorbital period as the 
X-ray radiation. 

Since the superorbital period shows considerable varia- 
tion (see, for example, Clarkson et al., 2003), folding upon 
the average period was not appropriate. The RXTE/ASM 
data were therefore folded using Table 2 in Trowbridge et al. 
(2007) which logs the times of superorbital zero-phase. This 
table continues up to the date of publication of their work, 
and so had to be extended up to make it possible to include 
all the optical data presented here. The RXTE/ASM X-ray 
data were smoothed using a five point moving average and 
the zero-phase points were estimated by eye. An updated 
edition of Fig. 5 from Trowbridge et al. (2007) showing the 
time evolution of the superorbital cycle length is presented 
here in Figure [7] 

The folded RXTE/ASM and OGLE light curves pro- 
duced as a result of folding the data at this varying super- 
orbital period are shown in Figure [8] The ASM light curve 
shows a clean shape (as would be expected from its L-S peri- 
odogram) with ~100% modulation. The OGLE light curve, 
on the other hand, is not so strongly modulated and only 
varies by ~3%. % 2 tests on the binned OGLE data were car- 
ried out: firstly using the ASM light curve as a model; then 
a straight line centered on the mean value. We calculated 
reduced x 2 values of 1.3 and 111.8 respectively, providing 
strong support for a correlated optical and X-ray superor- 
bital modulation. This is the first time such an effect has 
been noted. 
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Figure 8. Data folded at the super-orbital period - see text for 
details. X-ray data (upper panel) from RXTE/ASM in units of 
ASM counts/s. Optical data (lower panel) from OGLE-III I-band 
in flux units of 10 — 15 erg cm" 2 s _1 A~ 1 . 

5 DISCUSSION & MODELLING 

Looking first at the binary modulation of the I-band data 
(Figure Q we note the extra depth of the dip at phase 0.0 
(the time of the eclipse of the X-ray source) compared to 
phase 0.5 in this folded light curve. The difference in the 
depths of the dips is ~2-3%. Therefore one explanation for 
the discrepancy between the two minimima is that there is a 
small component of I-band light coming from the the accre- 
tion disk itself which is, of course, eclipsed at the same phase 
as the X-ray source. All the rest of the I-band modulation 
is due to the asymmetry in the shape of the supergiant star. 
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The fact that the visible band temperature is warmer dur- 
ing X-ray eclipse (see Figure [5]) adds further information, 
suggesting that the the source of the visible light in the ac- 
cretion disk is somewhat cooler than the surface of the BOI 
star - which is typically 30,000K. 

From the superorbital modulation shown in Figure [8] 
it is possible to determine that ~3% of the Fband light is 
modulated at the superorbital period. Since the nature of 
the superorbital modulation in the X-ray is attributed to 
the precession of a warped accretion disk (Clarkson et al., 
2003), it is therefore very likely that the source of visible 
light is also being modulated in the same manner. If this is 
the case, and the optical variation mirrors the 100% modu- 
lation seen in the X-ray profile, then one can make another 
estimate of the fraction of the light coming directly from the 
accretion disk, which is ~3%. Therefore both estimates of 
the fraction of Fband light arising from the accretion disk 
are in broad agreement at ~ 2 — 3%. Though the exact ori- 
gin of the optical superorbital modulation is still unclear, it 
worth noting that Rajoelimanana, Charles & Udalski (2011) 
suggest, based upon other OGLE data, that such long-term 
modulation may be present in many more SMC sources and 
hence may be more widespread than previously reported. 
In their paper they present an apparent correlation between 
orbital and superorbital periods in 17 systems (their Figure 
45) and suggest that the two periods could be linked by the 
size of the circumstellar disks around the Be stars and the 
light therefrom. Interestingly the two periods seen in SMC 
X-l do not place it comfortably upon an extension of their 
plot. This may be related to the suggestion that supergiant 
systems do not form circumstellar disks but only exhibit 
radial wind outflows. 

Rawls et al. (2010) recently redetermined the system 
parameters for six eclipsing high-mass X-ray binaries, in- 
cluding SMC X-l. The Eclipsing Light Curve (ELC) code 
of Orosz & Hauschildt (2000) was used to model all of the 
observational constraints in a self-consistent way. In the case 
of SMC X-l, the observational constraints include a U-band 
light curve from van Paradijs & Kuiper (1984), the semi- 
amplitude of the O-star's radial velocity curve K rmopt , the 
X-ray eclipse semiduration 6 e , the orbital period P, and the 
projected semimajor axis of the pulsar's orbit ax sin i (see 
van der Meer et al. 2007 and the cited references for a tabu- 
lation of these values) . Rawls et al. showed that the fit to the 
optical light curve was improved when an accretion disk was 
included in the model. In particular, the accretion partially 
eclipses the O-star, increasing the depth of the minimum 
near phase (the inferior conjunction of the pulsar) . 

We extended the analysis performed by Rawls et al. 
(2010) by including the mean phase-averaged OGLE light 
curve. When the model was extended to the /-band it pro- 
vided a very good fit to the mean OGLE light curve. An 
optimization using a Monte Carlo Markov Chain resulted in 
M ns = 1.145 + / - 0.074M©, M opt = 17.34 + / - O.9OM , 
and an inclination of i = 63.2 + / — 2.1°. Figure [9] shows 
the phased light curves and the best-fitting models. There 
are some small systematic differences near the quadrature 
phases in the fit to the OGLE light curve, possibly caused 
by some asymmetry in the accretion disk. In spite of these 
small deviations, the inclination of the binary is well con- 
strained, leading to well defined masses for the components. 

The neutron star mass is a little higher than that found 
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Figure 9. V & I-band data phase-folded at the orbital period 
and compared to synthetic lightcurves - see text for details. 



in Rawls et al. (2010), however, it confirms that it is very 
much on the low side for systems that have determined neu- 
tron star masses (see, for example, the sample in Rawls et 
al, 2011 and references therein). 
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